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Superconductivity at the interface between the insulators LaAlOs and SrTiOa has been tuned with 
the electric field efi^ect. The data provide evidence for a two dimensional quantum superconductor 
to insulator (2D-QSI) transition. Here we explore the compatibility of this phase transition line 
with Berezinskii-Kosterlitz-Thouless (BKT) behavior and a 2D-QSI transition. In an intermediate 
regime, limited by a finite size effect, we uncover remarkable consistency with BKT- criticality, weak 
localization in the insulating state and non-Drude behavior in the normal state. Our estimates for 
the critical exponents of the 2D-QSI-transition, z ~ 1 and 77 ~ 2/3, suggest that it belongs to the 
3D-xy universality class. 

PACS numbers: 74.78.-w, 74.40.-fk, 74.90.-fn, 74.78.Fk 



The conducting interface between LaAlOs and SrTiOa, 
two excellent band insulators, has been attracting a lot 
of attention[l, 2, 3, 4, 5]. Recently, two different ground 
states, magnetic Q and superconducting Q, have been 
experimentally identified. In a very recent report [8|, it 
was shown that the electric field effect can be used to map 
the phase diagram of this interface system. As the car- 
rier density is increased the system undergoes a 2D-QSI 
transition. A further increase reveals a superconducting 
dome. Moreover it was shown that the characteristics 
of the superconducting transition are consistent with a 
superconducting sheet of about 10 nm thick |9[. 

Here we attempt to unravel the nature of the phase 
transition line and its endpoint which separates the su- 
perconducting from the insulating ground state. For this 
purpose we explore the compatibility of the phase transi- 
tion with Berezinskii-Kosterlitz-Thouless (BKT) critical 
behavior [lOJ, [ll| and of its endpoint, where supercon- 
ductivity disappears, with a two dimensional quantum 
superconductor to insulator (2D-QSI) transition associ- 
ated with weak localization 12!|. In contrast to the disor- 
der tuned QSI transition, electrostatic tuning changes the 
carrier density without altering the disorder landscape 
[13| . In an intermediate temperature regime we uncover 
remarkable agreement with BKT-criticality in the super- 
conductor to normal state transition, weak localization 
in the insulating phase, and non-Drude behavior in the 
normal state. It is shown that both conduction mecha- 
nisms are limited by a finite size effect, whereupon the 
BKT-correlation length and the diverging length associ- 
ated with weak localization [12{ cannot grow beyond a 
limiting length, set by the extent of the homogeneous 
regions. Our analysis also reveals that the electrostatic 
tuned BKT-phase transition line ends at a 2D-QSI criti- 
cal point which appears to fall onto the universality class 
of the classical 3D-xy model. In addition we explore the 



Tc dependence of the vortex core radius and the vortex 
energy. These properties appear to be basic ingredients 
to understand the Tc variation. 

To explore the compatibility with BKT critical behav- 
ior we invoke the characteristic temperature dependence 
of the correlation length above Tc [ill , 



^ = Co exp (2^/ (6ii/2)) , t = \T/Tc - 1| 



(1) 



where ^o is the classical vortex core radius and b is 



14 



related to the energy needed to create a vortex 
ISl . Ild |. Note that b also enters the temperature de- 
pendence of the magnetic penetration depth A below 
the universal Nelson-Kosterlitz jump: A^ {Tc) /A^ (T) = 

1 -f 6 |i| /4) [lj|. Moreover, 5 is related to the vortex 



energy Ec in terms of 16l.ll7| 

b = 47TTyybn = f [Eel [kBTc)) . 



(2) 



Invoking dynamic scaling the resistance R scales in I? = 2 
as El 



Ro^C 



(3) 



where Zd is the dynamic critical exponent of the classical 
dynamics. Zd is usually not questioned to be anything 
but the value that describes simple diffusion: Zd = 2 [l9| . 
Combining these scaling forms and taking the occurrence 
of a 2D-QSI transition at Rqc = Ro {Tc = 0) into account 
we obtain 



with 



R{T)^Roexp{-bii{T-T, 



.1/2 



-1/2 



(4) 



bR = 4^T,^/ 76, ARo = Roc - Ro oc 1/^^, (5) 



because ARq is the singular part of Rq close to the 2D- 
QSI transition. To assess the compatibility with the char- 
acteristic BKT-behavior we analyze the data in terms of 



{d\nR/dTy^''^ = {2/hBf'^ (T - Te). 



(6) 



Fig. m shows {d\nR/dT) 



-2/3 



T for Vg = 40 V. Vg 
is the gate voltage allowing to tune the carrier density. 
In spite of the rounded transition there is an intermedi- 
ate regime revealing the characteristic BKT-behavior ([6|) . 
Noting that real systems are homogeneous over a limited 
domain only, the rounding may be attributable to a stan- 
dard finite size effect, whereupon the correlation length 
cannot grow beyond a limiting length, set by the lateral 
extent L of the homogenous regions [iSJ, |20|. In this 

case R{T) scales as R{T)exp(bR\T -Tc\~^^^] /Rq = 
f{x). fix) is the finite size scaling function with 
X =exp(bR\T-Tc\~^^^) /L2 cx £,^/L'^. li £, < L criti- 
cal behavior can be observed as long as f{x) ~ 1. This 
regime corresponds to the horizontal line in the inset of 
Fig. [TJ If ^ > L the scaling function approaches f{x) oc x 
so that R{T)/Rq tends to L~'^, the behavior indicated 
by the dashed line. In this context it should be recog- 
nized that the BKT correlation length does not exhibit 
the usual and relatively slow algebraic divergence as Tc is 
approached (Eq. (P)). For this reason the elimination of 
the finite size effect would require unprecedented sample 
homogeneity. The emerging agreement with BKT behav- 
ior, limited by a standard finite size effect, allows us to 
discriminate the rounded transition from other scenarios, 
including strong disorder which destroys the BKT behav- 
ior. It also provides the basis to estimate Tc, h^ (Tc) and 

— 1/2 

Rq (Tc), and with that b (Tc) and ^o oc ARq with rea- 
sonable accuracy. 

Applying this approach to the R (T) data for each 
gate voltage Vg we obtain the BKT-transition line de- 
picted in Fig. [21 displayed as Tc vs. i?n (T*), the nor- 
mal state sheet resistance at T* — 0.4K. We observe 
that it ends around R]jc — 4.28 kfi where the system 
is expected to undergo a 2D-QSI transition because Tc 
vanishes. With reduced i?n the transition temperature 
increases and reaches its maximum value, Tcm — 0.31 K, 
around R\j ~ 1.35 kil. With further reduced resistance 
Tc decreases. To identify deviations from Drude behav- 
ior [a oc n) in the normal state, which will be discussed 
later, we also included the gate voltage dependence of the 
normal state resistance. 

According to the scaling theory of quantum critical 
phenomena one expects that close to the 2D-QSI transi- 
tion Tc scales as 18ll21[ 



Tc 



(7) 



where S is the appropriate scaling argument, measuring 
the relative distance from criticality. V denotes the crit- 
ical exponent of the zero temperature correlation length 
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FIG. 1; (dlnR/dT)-^'-' vs. T for Vg = 40 V from 
Caviglia ei al. 8j where R = 3/5-Rn- The solid 
line is (dlni^/dT)"^''^ = 6.5 (T - Tc) yielding the esti- 
mates Tc = 0.27 K and (2/6^)^''^ = 6.5; the inset 
shows R exp(bR \T - Tcr'/^) /Ro vs. exp(bR \T - r,!"''"') 

with Ro — 1670 Q. The upper branch corresponds to 
T > Tc and the lower one to T < Tc. The solid line 
is {R/Ro)exp(bR\T -Tc\~'^^^) = 1 and the dashed one 

{R/Ro)exp(bR\T-Tc\-'/^^ = 0.4exp(fefl |T - TJ-^/^^ 

i^L)'. 
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FIG. 2; Tc vs. Ra (T*) (•)and Vg vs. Ra (T*) (•) at 
T* = 0.4 K from Caviglia et al. 8]. The sohd line is 
Tc = 1.17 ■ 10~''A7?(r* ) and the dashed one Vg = Vgc + 
1.39-10"^ Ai?3/2(r*) with Ai?(r*) ^{RaciT*)-RniT*)), 



Rdc (T*) = 4.28 kQ. and Vg, 



-140 V. 



^ (T = 0) oc S^'^ and z the dynanfic critical exponent. 
From Fig. [2] it is seen that the experimental data point 
to the relationship 



Tcoc Ai?(T*) oc AV;^/^ 



(8) 



close to quantum criticality. If the scaling argument is 



AR{T* 



1, while ii 6 = AVa 



2/3. Since 



the measured modulation of the charge density An2D in- 
duced by the gate voltage scales in the regime of interest 
as S 



/\V„ ex An2D oc T3/2, 



(9) 



we obtain zv ~ 2/3 if AVg or An2D are taken as scaling 
argument 5. To identify the correct scaling argument we 
use the fact that 5 oc AnsD holds ii{2 + z)V>2^. To 
check this inequality, given zF, we need an estimate of z. 
For this purpose we invoke relation ^, Roc — ^o (Tc) oc 
^(7^ (Tc) [zci = 2) which diverges as ^oiTc) oc ^ (T = 0) oc 
(5' 



oc Tf 



l/z 



so that the scaling relation 
Roc - Ro {Tc) oc ^0^2 ^Tc) oc (^2- 



rpljz 



(10) 

holds. Fig. [3] depicts the T^ dependence of fo C^c) « 
(i?Oc — -Ro (21;)) and 6, which is related to the vortex 

energy ii^c (Eq. (I2|)). Approaching the 2D-QSI-transition 
we observe that the data point to ^o (2^c) oc l/Tc, yield- 
ing for z the estimate z ~ 1 so that v ~ 2/3 with 
zv ~ 2/3. As these exponents satisfy the inequality 
(2 + z) 77 > 2 [221 we identified the correct scaling ar- 
gument, (5 oc An2D oc AVg. From the gate voltage de- 
pendence of the normal state sheet resistance shown in 
Fig. 2 it also follows that the normal state conductivity 
scales as Actq (T*) = ctq (T*) - a^c {T*) oc AVg^^. The 
2D-QSI transition is then characterized by the scaling 
relations 



Tc oc S'^/^ oc AR (T*) oc Ai?^^ (TJ oc AV^^/^ oc An2^ 
oc Aan (T*) ex Co^i (T,) , 



2/3 
2D 

(11) 



where Actq (T*) oc An^j^ reveals non-Drude behavior in 
the normal state. The product zv ~ 2/3 agrees with 
that found in the electric field effect tuned 2D-QSI tran- 
sition in amorphous ultrathin bismuth films |23j and the 
magnetic-field-induced 2D-QSI transition in Nbo.15Sio.85 
films [24] . On the contrary it differs from the value zF ~ 1 
that has been found in thin NdBa2Cu307 films using the 
electric-field-effect modulation of the transition temper- 
ature 25|. In any case our estimates, z ~ 1 and v ~ 2/3 
point to a 2D-QSI transition which belongs to the 3D-xy 
universality class [l^. Fig. [3] also depicts the T^ depen- 
dence of 6, which is related to the vortex energy. Since h 
tends to a constant in the limit T^ -^ 0, Eq. ^ implies 
db/dTc = and therewith Ec{Tc) oc ksTc, while the vor- 
tex core radius ^0 diverges as ^0 {Tc) oc l/Tc, in analogy 
to the behavior of superfluid "^He films [26|. The 2D-QSI 
transition is then also characterized by an infinite vortex 
core radius and vanishing vortex core energy. As Tc in- 
creases from the 2D-QSI transition, ^0 drops, while the 
vortex core energy increases. Finally, after passing the 
maximum Tc the vortex core radius continues to decrease 
with reduced Tc while b increases further. 

Having presented the evidence for an electric field effect 
tuned BKT line ending at a 2D-QSI-transition belonging 
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FIG. 3: Co {Tc) oc (7?oc - -Ro (Tc))"''"' (•) and h (•) vs . n 
from Caviglia et al. H where R = 3/57?n • The solid line is 
^0 {Tc) cc (.Roc - -Ro {Tc))~^^^ = 8 ■ 10-'' /Tc with Roc = 2700 

n. 
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FIG. 4; (jQ vs. T and (an (T) - guo) / {d\n{T)) vs. 1/ \ln{T)\ 



at Va 



-240 V from Caviglia et al. 



The solid line is 



o-Q (T) = o-Qo + dln(r) with aao = 1-2 • 10"'' {n~^) and d = 
1.2-10"^ (r2"^), the dashed one {an {T) - ano) / {din [T)) = 
1 and the dash-dot one (crn {T) — fgo) /d = 1.8 ex 1/L. 
The inset shows the magnetoconductivity ag vs. H, ap- 
plied perpendicular to the interface, at T = 0.03 K and 



K 



-300 V, from Caviglia et al. 



The solid line is 



an =4.51- 10"^ + 1.12 -10"^ In (JY) kO.'^ . 



to the 3D-xy universality class and non-Drude behav- 
ior in the normal state an important issue remains, the 
nature of the insulating phase. Fig. [3] shows the temper- 
ature dependence of ctq for Vg = —240 V, which is rather 
deep in the insulating phase [8| . In analogy to Fig. [1] 
we observe a rounded transition. In the present case 
the intermediate regime is compatible with weak local- 
ization, CTn {T) — CTno + rfln(r) [I4I. Given the evidence 
for the limiting length L, which prevented the attain- 
ment of the superconducting ground state, one suspects 
that it limits the growth of ^id oc i?n = I/ctdi the di- 



verging length associated with locahzation [l^, as well. 
In this case finite size scaling predicts that crn (T) scales 
as {(Jn{T) - (Tdc) I (c^ln ^)) = 9 iv) with y = S,id/L ex 
1/ {L [In (T)|). g{y) is the finite size scaling function and 
tends to 1 for y < 1. In this case the approach to the 
insulating ground state can be seen, while for y > 1 the 
crossover to g{y) — > y sets in and ctu^T) approaches the 
finite size dominated regime, where (Ju{T) — ctqo oc 1/L. 
A glance at Fig. |4] reveals that weak localization com- 
bined with the finite size effect describes the data very 
well. In addition, below Vg = —300 V d is indepen- 
dent of Vg and its value d = 1.2 ■ 10~^ Vt^^ is close to 
de = e^/ (tt/i) ~ 1.23 -10"^ r2"~\ generically attributed to 
electron-electron interaction [27| . The resulting evidence 
for weak localization is further substantiated by the ob- 
served negative magnetoresistance '8* and in particular 
by the high field behavior of the conductance depicted 
in the inset of Fig. [H There we observe consistency 
with the characteristic ln(i?) high field behavior [l2] and 
dau/d\n{H) ~ 1.12 • 10~^ri~^ is close to the theoreti- 
cal prediction dau/dXn [H) = e^ / (nh) ~ 1.23 • lO^^f^^i 
[28J . As a result, the failure to observe the superconduct- 
ing and insulating ground states directly, is attributable 
to a finite size effect, preventing the respective diverging 
lengths to grow beyond L. Nevertheless, the finite size 
scaling analysis and the high field magneto conductivity 
provide substantial evidence that these are the appropri- 
ate ground states in the homogenous and infinite system. 
In summary, we have shown that the electrostatically 
tuned phase transition line at the LaAlOa/SrTiOs inter- 
face, observed by Cavigiia et al. [8|, is consistent with a 
BKT-line ending at a 2D-QSI critical point with critical 
exponents z ~ 1 and V ~ 2/3, so the universality class 
of the transition appears to be that of the classical 3D- 
xy model. The normal state was shown to exhibit non- 
Drude behavior. To identify the superconducting and in- 
sulating ground states from the temperature dependence 
of the resistance we performed a finite size scaling anal- 
ysis because the growth of the finite temperature corre- 
lation length and the diverging length associated with 
localization turned out to be limited. Nevertheless, in 
the insulating state we observed in both, the tempera- 
ture and magnetic field dependence of the resistance, the 
characteristic weak localization behavior. Furthermore, 
we explored the Tc dependence of the vortex core radius 
and the vortex energy. These properties appear to be 
basic ingredients to unravel the nature of the variation 
of Tc- Approaching the 2D-QSI transition the vortex en- 
ergy tends to zero while the vortex core radius and the 
localization length diverge so the system is an insulator. 
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